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Abstract 

We analyze the uncertainties induced by different definitions of the momen- 
tum fraction z in the application of gluon to heavy quarkonium fragmentation 
function. We numerically calculate the initial g —>■ J ftp fragmentation func- 
tions by using the non-covariant definitions of z with finite gluon momentum 
and find that these fragmentation functions have strong dependence on the 
gluon momentum k. As \k\ —>■ oo, these fragmentation functions approach to 
the fragmentation function in the light-cone definition. Our numerical results 
show that large uncertainties remains while the non-covariant definitions of z 
are employed in the application of the fragmentation functions. We present 
for the first time the polarized gluon to J/ip fragmentation functions, which 
are fitted by the scheme exploited in this work. 



PACS number(s): 13.87.Fh, 14.70.Dj 



1 Introduction 



Fragmentation refers to the process of a parton which carries large transverse momen- 
tum and subsequently forms a jet containing the expected hadron pQ. At sufficiently 
large transverse momentum of the heavy quarkonium production, the direct leading order 
production scheme is normally suppressed while the fragmentation scheme becomes dom- 
inant, though it is formally of higher order in the strong coupling constant a s [2], [31 HI [5] . 

Generally, the fragmentation processes of heavy quarkonium H production can be 
expressed as [6] 

da[A + B -> H(p T ) + X] = V/ dx a [ dx b [ dzf a/A (x a ,Q) x 

~r Jo Jo Jo 

a,b,c 

fb/B(xb, Q)da(a + c + X)D c ^ H (z, Q), 

where a and b are incident partons in the colliding hadrons A and B respectively; f a /A 
and fb/B are the parton distribution functions at the scale Q 2 of the partonic subprocess 
a + b — > c + X; cis the fragmenting parton (either a gluon or a quark) and the sum runs 
over all possible parton contributions; D c _^h{z,Q) is the fragmentation function with 
respect to the scale Q 2 which can be obtained by evolving from the initial fragmentation 
function _D c _>#(z, Qq) using Altarelli-Parisi equations [7] 



Q^-D^ H (z,Q) = £ /' ^-P ij (z/y,Q)D j ^ H (y,Q) 
0(4 ■ Jz y 



where Pij are the splitting functions of a parton i into a parton j. The initial gluon 
fragmentation function D 9 ^h{z,Qo) is a universal function defined by factorization in 
the infinite momentum frame of the fragmenting gluon. It can also be obtained by cal- 
culating a specific physical process in perturbative QCD in the finite momentum frame 
of fragmenting gluon [U [S], where the z is defined in the Lorentz boost invariant form 
(e.g. Eq. ([I]) or The resulting initial fragmentation function, independent of the 

momentum of the parent gluon, is equivalent to the one derived in the infinite momentum 
frame. 

* = ^4- ( 2 ) 

In above equations we take the Z axis along the momentum of fragmenting gluon. Then, 
E H ,p*,E 9 and p 9 z are the energies and Z-components of the four-momenta of the frag- 
menting gluon and the produced heavy quarkonium H, respectively. In Eq. (pQ) z is defined 
as the usual light-cone form. 
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The Eqs. ([I]) and ([2]) are hard to be employed in the application of the gluon fragmen- 
tation functions, because they involve the transverse momentum of the resulting heavy 
quarkonium. Instead, usually the non-covariant definitions as follows are used approxi- 
mately: 



Ml + (p? ) 2 

5 W 



M^ + (pf) 2 +pf 

z = 4- (5) 
pi 

When the fragmenting gluon momentum \k\ — > oo, these non-covariant definitions are 
equivalent to the light-cone definition in Eq. (TjQ). However when the momentum of the 
fragmenting gluon remains finite, these non-covariant definitions of z may induce large 
uncertainties in the application of the fragmentation functions. The Eq. ([3]) can be re- 
expressed as 



|pf| = ^{zEoy-Ml, 

which shows that at certain values of z, there are two possibilities for momentum directions 
of quarkonium H. Similarly, the definition of Eq. (jlj) does not give unique direction of H 
for certain z, which does not comply with the original idea of fragmentation. Thus, the 
definition of Eq. ([5]) is better and widely used. About a decade ago, people noticed that 
the parton to heavy quarkonium fragmentation functions can be obtained analytically 
IH El El EE] ■ However, in fact the initial scale gluon fragmentation function can also be 
obtained numerically and works equally well as the analytic one in the application. 

In section 2, by using the perturbative QCD we numerically calculate the initial 
g — > J /if) fragmentation functions D g ^j/^(z, Mjm) with the non-covariant definitions 
given above at finite momentum of fragmenting gluon. We compare them with the frag- 
mentation functions in the light-cone definition Eq. ([1]). In section 3, We give out nu- 
merically the polarized gluon to J /if) fragmentation functions with light-cone definition 
of z, and transform it into analytic ones by fitting. We also present the x-square fitted 
functions for both polarized and unpolarized fragmentation functions in the light-cone 
definition of z. In section 4, we give a brief summary of our results. 



2 Unpolarized fragmentation functions w.r.t. differ- 
ent definitions of z and the resulting uncertainties 

The Feynman diagram of the process qg — * qg* — * J/ipgg is shown in Fig. [1] and the 
corresponding matrix element is composed of two parts: (1) the production of a virtual 
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Figure 1: The Feynman diagram of process qg — > qg* — > J j^gg 

gluon g* with momentum k and invariant mass y^; (2) the virtual gluon decaying into a 
J/ip and a gluon pair, i.e., 
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= -E^%(A)^M 2 >K(A), (6) 

where M( 1 £ A1 (A) and M2 (s)e*(A) represent the virtual gluon production and decay sectors, 
respectively. The A and k denote the virtual gluon and J/ip polarization indices. In 
deriving out the above second equation, the following condition on the summation of 
gluon polarization vector is used: 

v = -E £ m(aK(a). (7) 

A 

Then, the differential cross-section of process g + q — > q + g*\ g* — > J /if) + g + g, as shown 
in Fig. [H can be expressed as: 

da = \M\ 2 [d<f)} 

where the sum runs over the J /if) polarization for the unpolarized production, and 



da(qg^qg*) = S\ Pl + p 2 -k-p 6 )Y / \M x \l x [d^\ (2 ^ 9 , (9) 



dT K (g* ^ J/^gg) = |^**(* - Ps - Pa - p B ) ^ £ I^WIaa^]) , (10) 

= [d^] = -|^_ ) [dfc] = J] T^fer- (11) 

(27rJ- 1 2£; 6 i=3)4 5 (2vr) d 2^ 

Hence, the initial scale gluon fragmentation function corresponding to the finite gluon 
momentum satisfies 

^ = y y^ ^ ? y — ^ — • (i2) 
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Figure 2: The fragmentation functions w.r.t. Eq. (JTJ) (solid line) and ([2]) (dash line) 

Therefore, from Eq. (|T2|) one can extract the unpolarized initial gluon fragmentation 
functions with either boost-invariant or non-invariant definitions of z. We numerically 
calculate the fragmentation functions with various definitions of z with a specific gluon 
momentum and compare the differences among them. 

Fig. [2] is the fragmentation functions corresponding to the boost-invariant definitions 
of z given in Eq. ([1]) and (j2J), which are numerically calculated using the FDC (Feynman 
Diagram Calculation) program [11] and they agree with those given in Refs. [I] and 
[12] , respectively. Given different virtual gluon momentums in the numerical calculation, 
we obtain the same fragmentation function distribution versus z as shown in Fig. El 
which indicates explicitly the independence of the fragmentation functions on virtual 
gluon momentum, as they should be from the definitions of z. 

Figs. [3] and H] present the fragmentation function distributions over variable z in 
non-covariant definition of Eqs. ([3]) and (j4l). It shows that the fragmentation function 
distributions rely on the momentum \k\ of the virtual gluon. When \k\ — > oo, the cor- 
responding fragmentation functions approach to the one in the light-cone definition of 
Eq. 01) . However, when | A: J is finite, these fragmentation functions distribute distinctively 
from the one in the light-cone definition, which will bring certain uncertainties in the 
application of the fragmentation functions. With the non-covariant definition in Eqs. (J3]) 
and (jlj), there is the possibility that the J/ip momentum is opposite to the virtual gluon 
momentum k, and this possibility becomes larger as \k\ decreases. Therefore, the large 
peaks at small z in Figs. ([3]) and (j4j) correspond to the lower momentum \k\. 

Figure [5] gives the distribution of the fragmentation functions versus z in the definition 
of ([S]) at various virtual gluon momentums. Since it is possible that the J/ ijj may possess 
momentum opposite to that of the virtual gluon, in the definition of ([S]), the condition of 
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Figure 3: Fragmentation functions w.r.t. Eq. 




Figure 5: Fragmentation function w.r.t. Eq. (jSJ) 



z < exists. And because of the non-zero invariant mass of the virtual gluon, even z > 1 
happens. Thus, in Fig. [5] the fragmentation functions distribute in the scope of — 1 < z < 
2. As shown in the figure, when the virtual gluon momentum \k\ —>■ oo, fragmentation 
functions approach to the one in the light-cone definition. When \k\ = 5, 10,20,80 and 
1280GeV, the ratios of fragmentation probabilities in the definition of (5), that is the 
integrations of the corresponding fragmentation functions over z, to the one in light-cone 
definition are 0.49, 0.65, 0.77, 0.91 and 0.97 respectively in the range of < z < 1. This 
indicates that the smaller the virtual gluon momentum is, the more possible the produced 
J/ip will be opposite to the gluon. Therefore, uncertainties occur in the application of 
the fragmentation function in the definition of ([5]) when gluon momentum \k\ is not big 
enough. While the virtual gluon energy is greater than 80GeV, the uncertainties induced 
by the definition of §5§ will drop to less than 10%. 

Figure [6] presents the differential cross-sections of processes p p(p) —>■ g g versus gluon 
energy at the colliding energies of Fermilab Tevatron and LHC respectively. In our nu- 
merical calculation, the parton distribution function of CTEQ6L [13] is employed. As 
shown in the figure, the cross-section decreases quite fast as the gluon momentum in- 
creases. Whereas, from the analysis above, we know that the at some low energies of 
the fragmenting gluon there are large discrepancies between the fragmentation functions 
with non-covariant definitions and boost-invariant definitions of z. Therefore, it is possi- 
ble that the non-covariant definitions of z may induce large errors in the application of 
fragmentation functions. 
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Figure 6: The energy dependence of the differential cross-sections of p p(p) — ^ g g process 
at LHC and Tevatron. The horizontal variable E denotes the final gluon energy. 

3 Polarized fragmentation functions w.r.t. light-cone 
definition of z 

The polarized fragmentation functions of gluon fragmenting into P-wave charmonium 
Xcj exist in the literature [14J, however, there has been no relevant work on polarized 
fragmentation functions of gluon to S-wave vector Quarkonium states. With the FDC 
program and procedure described in above, in fact it is straightforward for us to obtain 
the gluon to polarized J/i/j fragmentation functions numerically. And then get the analytic 
ones by fitting, which are equivalent in use as the ones obtained directly from the analytic 
calculation, at least to large extent. 

With the numerical method described in Sec. 2, we can also obtain the polarized 
fragmentation functions which satisfy 



'^UMf«) = i*»r J^i ^"2 m "\ <i4) 



o '-* V ' ;J '" - J Ju} /t wsW J T tb 

where D^j,^ and D^j^ are fragmentation functions of gluon to J /if) in transverse 
and longitudinal polarizations, respectively. The polarized fragmentation functions can 
then be extracted from Eqs. (1131) and (j!4p in different definitions of z. Fig. [7] shows the 
polarized fragmentation functions corresponding to the light-cone definition of Eq. (TO . 
Explicit numerical calculation indicates that these polarized fragmentation functions in 
light-cone definition are independent of the virtual gluon momentum, as they should be. 
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In Ref. [8] , the initial unpolarized fragmentation function of gluon to J / ip is given in a 
form of two-dimensional integrals, which in practice can only be evaluated numerically and 
thus is not very convenient for the application of the Altarelli-Parisi evolution equation 
[7j. In use of the obtained initial gluon to unpolarized J /if) fragmentation function [8], 
one can integrate out the integrals numerically p~5] and fit the z distribution into a 
combination of functions of z. Then the application of the gluon to J/if> fragmentation 
function will be more easier for the phenomenological use. To our knowledge the polarized 
fragmentation functions of gluon to J/ if) are still absent so far in literatures. However, 
they may be necessary for future careful analysis about the polarization situation of 
Quarkonium production. 

We proceed the x _sc l uare fitting of gluon to polarized and unpolarized J /if} fragmen- 
tation functions into polynomials and exponents in the light-cone definition of z, which 
tells 

= 1 e^-'\l-z) E On*", (15) 

n=l,9 

\2 m J/i>) n=l,9 

Here, a s (Mj/^), Mj/$ and -R(O) are strong coupling, the J/if> mass and radial wave 
function at the origin, respectively. The fitted coefficients of a n and a„ are presented in 
Tabled 



a 


-17.0402262 


T 
% 


-17.3527537 


ax 


5.1820321 x 10" 7 


T 

a{ 


3.3809309 x 10~ 7 


a 2 


-7.5744212 x 10~ 6 


T 


-5.0184819 x 10~ 6 


a 3 


1.2310370 x 10" 4 


T 

aj 


8.2006550 x 10" 5 


a 4 


-7.6536059 x 10~ 4 


T 

a\ 


-5.1990241 x 10" 4 


a 5 


3.1398595 x 10~ 3 


T 


2.1421231 x lO" 3 


a 6 


-7.3523074 x 10~ 3 


T 

a 6 


-4.9941798 x 10~ 3 


a-j 


9.2419267 x 10" 3 


T 


6.2434863 x 10~ 3 


a 8 


-5.8527342 x 10~ 3 


T 

a£ 


-3.9354343 x 10" 3 


a 9 


1.4728577 x 10~ 3 


T 
«9 


9.8679226 x 10~ 4 



Table 1: Coefficients of fitted functions for the unpolarized and transversely polarized 
initial gluon to J /if) fragmentation functions 
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Figure 7: Polarized fragmentation functions w.r.t. to Eq. [3] 

4 Summary 

In this work, we show that although the fragmentation functions corresponding to the 
boost invariant definitions of (DO) and (j2j) are equivalent to the ones in the infinite mo- 
mentum frame of gluon and thus independent of fragmenting gluon momentum, they are 
inconvenient in use, because they involve the transverse momentum of the parent gluon. 
Instead, the non-covariant definitions, such as in Eqs. d3J), (jlj) and (J5j), are used as an 
approximation. They are equivalent to the definition of the light-cone coordinate form 
when the fragmenting gluon momentum \k\ — ► oo. However, as shown in our calculation, 
these non-covariant definitions may induce large errors while \k\ is finite, especially small. 
In practice, the definition of ([3]) and (jl]) are somewhat inconvenient because they possess, 
in the range of < z < 1, the possibility of J/ip moving in the opposite direction of the 
fragmenting gluon. The definition of ([5]) excludes this situation, but it may bring larger 
errors in the colliding energies of Tevatron and LHC This should be taken into account 
in the careful phenomenological analysis in the application of the gluon fragmentation 
functions. 

In this paper we give out the initial fragmentation functions of gluon into the polarized 
J I ip for the first time. We use a different method in getting them from the normal analytic 
calculation in the literature. The x _sc L uare fittings for both polarized and unpolarized 
fragmentation functions of gluon into J/ip are done and presented, which are suitable 
for future phenomenological use. Finally, although the polarized fragmentation functions 
obtained in this work are schematically for charmonium, the J/ip, in fact they can be 
directly applied to the T system with some simple replacements. 
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